Bures metric over thermal state manifolds and quantum criticality 
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We analyze the Bures metric over the manifold of thermal density matrices for systems featuring 
a zero temperature quantum phase transition. We show that the quantum critical region can be 
characterized in terms of the temperature scaling behavior of the metric tensor itself. Furthermore, 
the analysis of the metric tensor when both temperature and an external field are varied, allows 
to complement the understanding of the phase diagram including cross-over regions which are not 
characterized by any singular behavior. These results provide a further extension of the scope of 
the metric approach to quantum criticality. 
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I. INTRODUCTION 

These years are witnessing an increasing research ef- 
fort at the intersection of Quantum Information Sci- 
ence [l| and more established fields like theoretical con- 
densed matter Q. It belongs to this class the approach 
to quantum phase transitions (QPT) (3) based on the 
information-geometry of quantum states that has been 
recently proposed in Ref Q and Further develop- 
ments, for specific, yet important, class of quantum states 
have been then reported in 0, i, ©, M, M, [H [H [3 . 
The underlying idea is deceptively simple: the major 
structural change in the ground state (GS) properties 
at the QPT should reveal itself by some sort of singu- 
lar behavior in the distance function between the GSs 
corresponding to slightly different values of the coupling 
constants. This intuition can me made more quantitative 
by analyzing the leading order terms in the expansion of 
the quantum fidelity between close GSs. 

A general differential-geometric framework encompass- 
ing all of these result has been offered in Ref [l5| ■ There 
it has been shown that these leading order terms do cor- 
respond to a Riemannian metric g over the parameter 
manifold. This metric g is nothing but the pull-back of 
the natural metric over the projective Hilbert space via 
the map associating the Hamiltonian parameters with 
the corresponding GS. In the thermodynamical limit the 
singularities of g correspond to QPTs. In Ref. [lj| the 
nature of this correspondence has been further investi- 
gated and it has been shown that both the metric ap- 
proach to QPT and the one based on geometrical phases 
[171 . [HI ] can be understood in terms of the critical scaling 
behavior of the quantum geometric tensor [T^ |. 

The conceptually appealing and potentially practically 
relevant feature of this strategy consists of the fact that 
its viability does not rely on any a priori knowledge of the 
physics of the model e.g., order parameters, symmetry 
breaking patterns,... but just on a universal geometrical 
structure (basically the Hilbert scalar product). Very 
much in the spirit of Quantum Information the metric 
approach is fully based on quantum states rather than 



Hamiltonians (that might be even unknown), once these 
are given the machinery can be applied. 

In this paper we further extend the scope of this metric 
approach by considering the manifold of thermal states 
of a family of Hamiltonians featuring a zero-temperature 
PT. In [20] it was shown that by studying the mixed-state 
fidelity [21| between Gibbs states associated with slightly 
different Hamiltonians one could detect the influence of 
the zero-temperature quantum criticality over a finite 
range of temperatures. Here we will refine that anal- 
ysis and make it more quantitative by resorting to the 
concept of Bures metric between mixed quantum states. 
This metric provides the natural finite-temperature ex- 
tension of the metric tensor g studied in the GS case and 
corresponds again to the leading order in the expansion of 
the (mixed-state) fidelity between close states i.e., associ- 
ated with infinitesimally close parameters. By analyzing 
the case of the Quantum Ising model we shall show how 
the quantum-critical region above the zero-temperature 
QPT can be remarkably characterized in terms of the 
scaling behavior of the Bures metric tensor. 

The paper is organized as follows: in sect. II we intro- 
duce the basic concepts about mixed-state metrics and in 
Sect III we specialize them to the case of thermal (Gibbs) 
states. In sect. IV we provide generalities about quasi- 
free fcrmion systems and in Sect VI we analyze in detail 
the Bures metric tensor for the quantum Ising model. 
Finally in sect V conclusions and outlook are given. 



II. PRELIMINARIES 

The Bures distance between two mixed-states p and a 
is given in terms of the Uhlmann fidelity [2l| 



F(p, a) = trvV/Vp 1 / 2 



(1) 



by d B (p, a) = y/2[l-F(p,(T)]. 

The starting point of our analysis is provided by the 
following expression for the Bures distance between two 
infinitesimally close density matrices (see e.g. [22] for a 



2 



derivation) 



ds 2 (dp) := d% (p, P + dp) = ^^2 



\{m\dp\n)Y 



P„ 



Pn 



(2) 



where \n) is the eigenbasis of p with eigenvalues p n i.e., 
p = J2nPn\ n ){ n \- Even though in the sum in ((2|) p n and 
p m cannot be simultaneously in the kernel of p, since 
|n), |m) S Ker(p) (n|c?p|m) = 0, one can formally ex- 
tend the sum to all possible pairs by setting to zero the 
unwanted terms. For p pure i.e., p — one has 

dp = \dip)(ip\ + \ip)(dip\ from which one sees that the di- 
agonal matrix elements of dp are vanishing and one is left 

With ds% =E meK er(p) IW|m)| 2 = (#|(l-h^|)|#). 

This expression coincides with the Ricmannian metric 
considered in The bures metric (J2]) is tightly con- 

nected to the so-called quantum Fisher information and it 
appears in t he q uantum version of the celebrated Cramer- 
Rao bound [23| . This suggests the possible relevance of 
the results that we are going to present in this paper to 
the field of quantum estimation [24| . 

To begin with we would like to cast Eq. ^ in a from 
suitable for future elaborations. Let us first differenti- 
ate the density matrix dp = J2n(dPn\n)(n\ +p n \dn){n\ + 
p n \n)(dn\) and consider to begin the matrix element 



(dp)ij. We observe that = 5ij - 
whence (i\dp\j) = Si jdpi + (i\dj)(pj 
expression back into ([2|) one obtains 



(di\j) = -ii\dj); 
Pi). Putting this 

= - V ^ + - J2 \(n\dm)\ 2 ^" ~ P ™ )2 . (3) 



4^ p n 



n^m 



Pn + Pm 



This relation is quite interesting since it tells apart 
the classical and the quantum contributions. Indeed the 
first term in ([3]) is nothing but the Fisher-Rao distance 
between the probability distributions {p n } n and {p n + 
dp n }n whereas the second term takes into account the 
generic non-commutativity of p and p' := p + dp. We 
will refer to these two terms as the classical and non- 
classical one respectively. When [p 1 , p] = the problem 
gets effectively classical and the Bures metric collapses 
to the Fisher-Rao one; this latter being in general just a 
lower bound (HHl- 

Before moving to the analysis of the metric @ we 
would like to comment about the connection with the 
recently established quantum Chernoff bound [26] . This 
latter, denoted by £,qcb, is the quantum analogue of the 
Chernoff bound in classical information theory; it quan- 
tifies the rate of exponential decay of the probability of 
error in discriminating two quantum states p and a when 
a large number n of them is provided and collective mea- 
surements are allowed i.e., P err ~ exp(— u^qcb)- The 
Chernoff bound naturally induces a distance function 
over the manifold of quantum states with a well defined 
operational meaning (the bigger the distance between 
the states the smaller the asymptotic error probability 
in telling one from the other). In 26] it has been proven 
that exp(-^Qcs) = mm < s <i tr (p s cr 1_s ) < T(p,a) and 



that for infinitesimally close states i.e., a = p + dp, one 
has 



d« 2 

as QCB 



\{m\dp\n)Y 



(Vp^+Vp^) 2 ' (4) 



From this expression we see that the distinguishabil- 
ity metric associated with the quantum Chernoff bound 
has the same form of the Bures one (J2]) but the de- 
nominators p n + p m are replaced by {^/pn~i + ^/Pn) 2 - 
Using the inequalities {^/PnZ + y/p^) 2 > Pn + Pm and 
2(Pn + Pm) > {\fp~m + \fp~n) 2 one immediately sees that 



~Y < ds 2 QCB < ds 2 . 



(5) 



This relation shows that, as far as divergent behavior 
is concerned, the Bures and the Chernoff bound metric 
are equivalent i.e., one metric diverges iff the other does. 
On the other hand in the metric approach to QPTs the 
identification of divergences of the rescaled metric tensor 
and their study plays the central role [l5(. Therefore 
one expects the two distinguishability measures to convey 
equivalent information about the location of the QPTs. 
Though most of the calculations that are reported in this 
paper could be easily extended to the Chernoff bound 
metric, here we will limit ourselves to the analysis of the 
Bures metric ©. 



III. THERMAL STATES 

From now on we specialize our analysis to the case of 
thermal states. If the Hamiltonian smoothly depends on 
a set of parameters, denoted by A, living in same man- 
ifold M. one has the smooth map (A, /3) — > p{fl, A) := 
z -i e -f)H(\)^ ^ Z _ tve -l3Hy wh at we are going to study 

in this paper is basically the pull-back onto the (A, (3) 
plane of the Bures metric through this map. This is the 
obvious finite-temperature extension of the ground-state 
approach of Ref . [ijl . 

We start by studying the Bures distance when T ^ is 
fixed and for infinitesimal variations of the Hamiltonian's 



parameters A. Notice first that p = Z 1 



-0E n 



n)(n\ 



where E n and |n) are the eigenvalues and eigenvec- 
tors of the Hamiltonian operator H. With a standard 
reasoning, by differentiating the Hamiltonian eigenvalue 
equation one finds that (i\dj) = (i\dH\j)/(Ei — Ej). 
Moreover one easily sees that dpi = d{e~@ Ei / Z) = 
—/3pi(dEi — (J2j dEjPj)), therefore the first term in equa- 
tion ^ can be written as /3 2 /4 J2 t P%{dE 2 — (dE)) 2 where 
{dE)p := J2j dEjPj). This means that the Fisher-Rao 
distance is expressed by the thermal variance of the di- 
agonal observable dHd := • dEj\j)(j\ times the square 
of the inverse temperature. Summarizing 



ds% 



{{dH 2 ) p - (dH d ) 2 ) + 



3 



2 ^ 



(n\dH 



TO 



E n — E„ 



2 (e-^~ - < 
Z(e-^» 



(6) 



The two terms correspond to the first and second term 
of ([3]) respectively and they depend on (3 and on the 
other parameters of the Hamiltonian. For example, when 
a single parameter h is considered, the Bures distance 
defines a simple metric that can be expressed in term of 
the classical and non-classical part 



(7) 



such that ds 2 B — ghh(h, (3)dh 2 . Let us now explore the 
behavior of the Bures distance in presence of infinitesimal 
variations of both the temperature (/3 variations) and a 
field h in the Hamiltonian. It is easy to see that the 
variation of (3 only affects the Fisher Rao classical term 
in ©. In fact the variation dH in ©, or analogously the 
variations \dm) in ([3]), are taken with respect to h only. 
The calculations can be summarized as follows. We first 
have to expand the dp n as dp n = (d/3p n )dp + (dhp n )dh. 
We have that 



and 



(d(3 Pn )df3 = p n [(H) - E n ] df] 



(d h p n )dh = @p n [(d h H d ) - d h E n ] dh 



where E n — E n (h). The complete classical term of 
the Bures distance can be written expanding (dp„) 2 — 
(dhPndh) 2 + (d/3p n d(3) 2 + 2d0p n d h p n df3dh, and summing 
over n. We thus have three different contributions: 



[dp r , 



\ { [(H 2 ) (H) 2 ] dp 2 



dh 2 



+p 2 [{[(d h H) d } 2 ) - ((d h H) d ) 2 

+2/3 [(H(d h H) d ) - ((d h H) d )(H)} d(3dh} 

(8) 



These terms correspond to the elements of the metric 
9/3,/3) 9h,/3i 9h h respectively. The full metric can be 
written once one calculates the non classical term in ([3]) . 
The infinitesimal Bures distance can then be written in 
terms of the 2x2 metric tensor g as : 



ds 2 



{dh, dP) g 



dh 

dp 



fJ 



9hh gh0 

ghp g/3/3 



(9) 



where again g hh {h, P) = g c hh (h, 0) + gffi(h, P). 

It is at this point interesting to check whether, for 
p — ► oo , one recovers the known results for ground- 
state (pure) fidelity and metric tensor. In order to 
do that we will consider separately the classical and 
non-classical term in ([3]). In fact \\p(P) — p(oo)\\i = 



(i - po) + E 



n>0 



Pr, 



< 



2£ 



n>0 ' 



-P(E n -E ). f rom w hich 

one sees that, for finite-dimensional systems, the ther- 
mal density matrix converges (in trace norm) exponen- 
tially fast to the projector over the ground state |0). 
Then it follows that all the expectations values will con- 
verge exponentially fast to their zero-temperature lim- 
its: |tr(Ap(/3))-tr(A P M)| < \\A\\\\p(P) - pMlli, this 
in turn guarantees that the covariances of diagonal op- 
erators in the Fisher-Rao term ([8]) are vanishing (since 
e.g., dH d is diagonal (dH 2 )^ = (dH) 2 ^.) in the zero- 
temperature limit. In the infinite dimensional case the 
convergence to zero of this term will typically be only 
algebraic in the region where the smallest excitation gap 
is small compared to the temperature, whereas it will be 
exponential elsewhere. The overall convergence behav- 
ior of the classical term for /3 — > oo depends now on the 
detailed interplay between the decay of covariances we 
just discussed in (jHJ) and the divergence of the powers of 
P in front of them. An analysis of the zero-temperature 
limit of these terms will be provided later for the quan- 
tum Ising model. We will see that all the classical terms 
vanish in the zero temperature limit but at the critical 
value of the parameter. As far as the second non-classical 
term in ([3]) is concerned one has just to notice that from 
lim^oo p n (P) — S n .o it follows that the only contribu- 
tions will come from the elements involving the ground- 
state i.e., (0\dj) = {0\dH\j)/(Ej - E a ). This completes 
the remark. 

Before moving to the next sections, where we will spe- 
cialize the previous results to the particular case of the 
Quantum Ising model, we would like to notice that the 
variation of the Bures distance with temperature only, 
given by the element gpp of the metric, is precisely pro- 
portional to the specific heat c v [l5j . i.e. 



ds 2 B = ^((H 2 ),- 



dp 2 



This simple fact was already observed in [l5[ and 
and provides, we believe, a neat connection between 
quantum- information theoretic concept, geometry and 
thermodynamics. 



IV. QUASI FREE FERMIONS 

In this section we specialize the study of the behav- 
ior of the Bures metric to systems of quasi-free fermions 
when one has the variation of one parameter h of the 
Hamiltonian and of the temperature T. The results that 
we present here are a finite-temperature generalization of 
those given in Refs Q and [|[ and directly related to the 
mixed-state fidelity ones reported in [20 ]. 

The quasi-free Hamiltonians we consider are given, af- 
ter performing a suitable Bogoliubov transformation, by 



H 



(10) 



4 



where A„ > and r\ v denote the quasi-particle energies 
and annihilation operator respectively. One has that v 
is a suitable quasi-particle label, that for translation- 
ally invariant systems amounts to a linear momentum; 
the ground state is the vacuum of the r\ v operators i.e., 
X] V \GS) = 0, Vv. The dependence on the parameter h is 
both through the A^'s and the rj v 's. 

We now derive the explicit general form of the Bures 
distance ([2]) starting from the classical part (8J. We ob- 
serve that the (many-body) Hamiltonian eigenvalues are 
given by Ej = ^2i V n v h. v where the n„'s are fermion oc- 
cupation numbers i.e., n v = 0, 1. Therefore we have that 
dEj = ^2 v n v dK v and (dEj)p = J2 U { n v)pdA v where the 
averages are easy to compute since the probability distri- 
bution of the dEj factorizes over the z/'s. Furthermore, 
{n^riv) p- (rip) p(n u ) p = 6^, v (n v )p (1 - {n v )p) and we can 
thus write: 



(d Pn ) 2 



I^(n fc )(l-(n fc ))x 

{A 2 d(3 2 

+p 2 {d h A k f dh 2 
+2/3A k d h A k dpdh}. (11) 



The term in dh 2 is the classical term due to the in- 
finitesimal variations of the parameters of the Hamil- 
tonian at fixed T and it corresponds to the variance, 
see ©, var(Hd) = Y.A^ffi 1 ~ {n v )p)dAl. Since we 
are dealing with independent free-fermions one has that 
(n v )/3 = (exp(/3A 1/ ) + l) -1 , whence 



dsi = '— 



(d h A k f 



16 cosh 2 0^/2) 



dh 2 



(12) 



In order to compute the non-classical part of Eq. ([3]), 
one has to explicitly consider the eigenvectors of (fTT))) . 
Following the notation of Ref. [7| one has \m = 

{a,,,a_„]»o) = ®v>o\a v ,a- u ) where 

|0„0„) = cos(^/2)|00),,_, y -sin(^/2)|ll),,_, y , 
\0 V 1- V ) - |01>,,_„ |1„0_,,) - \10) v ,- v , 
\l v l„) = cos(f„/2)|ll) I/) _ 1/ + sm(^/2)|00) I/ ,_ 1/ . 

We assume now that parameter dependence is only in 
the angles #„'s (this assumption holds true for all the 
translationally invariant systems). It is easy to see from 
the above factorized form that the only non vanishing 
matrix elements (n\dm) are given by (0 u 0- u \d\l u l^ u ) = 
d9 v /2 and that the thermal factor (p n — p m ) 2 /{Pn + Pm) 
has the form sinh 2 (/3A l/ )/[(cosh(/3A ly ) + l)(cosh(/9A I/ )] = 
(cosh(/3A ly ) — l)/cosh(/3A„). Putting all together one 
finds 



ds r 



1 ^ cosh(/3A„) - 1 
4 ^ cosh(/3A„) 



(d h e v fdh 2 (13) 



u>0 



We finally note that the two elements (fT2|) and (fl~3"l) define 
the metric element ([7]). The results of this section can be 
applied to any quasi- free fermionic model (|10|) . 



V. QUANTUM ISING MODEL 

We are now going to discuss in some detail the behavior 
of the metric tensor for a paradigmatic example in the 
class of quasi-free fermionic models, the Ising model in 
transverse field. The model is defined by the Hamiltonian 



■hot 



(14) 



At T = this system undergoes a quantum phase tran- 
sition for h = 1 . For h < 1 the system is in an or- 
dered phase as the correlator (af <t*)t=o tends to a non 

zero value: lim r _ KX> (of of )t=o = (l — h 2 ) 1 ^ 4 . The exci- 
tations in this region are domain walls in the a x direc- 
tion. Instead for h > 1 the magnetic field dominates, 
and excitations are given by spin flip over a paramag- 
netic ground state. The transition point h = 1 is de- 
scribed by a c = 1/2 conformal field theory, which im- 
plies that means that the dynamical exponent 2=1; 
the correlation function exponent is v = 1. As is well 
known :[3J , a signature of the ground state phase diagram 
remains at positive temperature. In the quasi classical 
region T -C A, where A = |1 — h\ is the lowest excita- 
tion gap, the system can be described by a diluted gas 
of thermally excited quasi-particles, even if the nature of 
the quasi-particles is different at the different sides of the 
transition. Instead in the quantum critical region T> A 
the mean inter-particle distance becomes of the order of 
the quasi-particle de Broglie wavelength and thus quan- 
tum critical effects dominate and no semiclassical theory 
is available. In each of the above described regions of the 
(h,T) plane the system displays very different dynamical 
as well thermodynamical properties. For example, in the 
quantum critical region the specific heat approaches zero 
linearly with the temperature (this is in fact a general 
feature of all conformal field theories), whereas in the 
quasi-classical regions the approach is exponentially fast. 



A. Bures metric tensor in the (h,T) plane 

We now investigate whether the signature of the physi- 
cally different regions can be revealed by analyzing the el- 
ements of the metric tensor defined by the Bures distance. 
We begin by studying the temperature dependence of the 
metric tensor when only the external field is varied i.e., 
the term g^h (h,T), see Eq. Q. The Hamiltonian CH]) 
is equivalent to a quasi-free fermionic model, and fol- 
lowing our previous notation one has e k — cos (fc) — h, 
A fe = sin(fc), A k = y/e 2 + A 2 and tan(tf fc ) = A k /e k . 
Using formulae (p~2|) and (fT3|) it is straightforward to write 
([7]). After rescaling g — > g/L and passing to the thermo- 
dynamic limit we obtain 



9hh 



I6n 



1 



ff cosh(/3A fc ) 



dk 
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nc 

9hh — 



Stt 



cosh(/3A fc ) - 1 A? 
cosh(/3A fc ) A* d 



first given in Q which we re-write here as an energy 
integral 



The integrals are better evaluated by transforming mo- 
mentum integration to energy integration in a standard 
way. As previously noticed, on general grounds, the clas- 
sical term vanishes when the temperature goes to 
zero. In the quantum-critical region f3A rs 0, and one 
obtains the following low temperature expansion: 



9hh 



96h 2 



T + O (T 2 ) 



Instead in the quasi-classical region where j3A ^> 1 the 
fall-off to zero is exponential. With a saddle point ap- 
proximation one obtains 



9hh 



32nh 



T -3/2 e -A/T + lower order ^ 



We now analyze the scaling behavior of the non- 
classical term of the metric g. From the results of 0, 0] 
it is known that the geometric tensor at zero tempera- 
ture diverges as A -1 when A — > 0. The non classical 
term matches this ground state behavior from positive 
temperature. Indeed, in the quantum-critical region the 
integral is well approximated by 



nc 
9hh 



8irh 2 



2 " cosh(x)-l\/W 2 
cosh (x) x 



-dx. 



For large (3 (low temperatures) this expression can be 
Laurent expanded and the resulting integrals can be re- 
summed using residuum theorem, giving 



n nc — — 

9hh - h 2 



^T- l -±- + 0{T) 

TT Z 16 



(15) 



where C is Catalan's constant C — 0.915966 . . .. 

We would like to point out that the behavior of the 
metric tensor in the quasi-critical region, can be inferred 
from dimensional scaling analysis in much the same spirit 
as was done in 16] for the zero temperature metric ten- 
sor. From Eq. (O we see that the scaling dimension of 
g^h is A nc — 2Ay — 2z — d, where Ay is the scaling di- 
mension of the operator dH, z is the dynamical exponent 
and d is the spatial dimensionality. Following [l6| (/? now 
plays the role of the length) we obtain 



nc 
9hh 



T 



(16) 



In the present case, z = d = Ay = 1 (the scaling dimen- 
sion of of -a free fermionic field- is one) which agrees 
with p5|) . 

We now pass to analyze the behavior of g^ in the 
quasi-classical region i.e. when 0A ^> 1. In this case the 
"temperature" part of the integral is never effective, i.e. 
one has 1 ~ 1, so it is quite clear that, in first 

approximation, one recovers the zero temperature result 



nc 
9hh 



8irh 2 



\i+h\ 
\i-h\ 



(h-iy 



(i + hf 



dus, 
(17) 

and we assumed h > 0. For small values of the gap, - 
hence we are in a situation where we consider first the 
limit T — > and then A — > — we observe the following 
divergence 

ff j£(T = 0,A-0) 



16A' 



which is a result also reported in [1, 0|- Instead 
when the gap is large, — so that we are necessary on 
the h > 1 side — we can approximate the radical in 
Eq. (fl7| with an ellipse centered at (h, 0) with semi- 
axes r. r = 1 and r„ — 2h, that amounts to write 



(h - I) 1 - cj 5 



-(i + hy 



2hJl-(co-hy 



In this case the integral gives 
<7^(T = 0,A»1) 



8^/2^ _l) 3 /2 8A4- 

Again, by doing a saddle point approximation one re- 
alizes that the zero temperature results are approached 
exponentially fast with the temperature, more precisely 
one has 



■3/2 -A/T 



gll G9A » 1) = gll (T = 0) — const, x T 



We now extend our analysis to the other terms of met- 
ric tensor ©. When we consider the case in which both 
the temperature and the field h are varied two new ma- 
trix elements come into play: 



9tt 



9hT 



16tt J_„ cosh(/3A fe ) + 1 
167ri_ w cosh(/3A fe ) + l 



-dk 



dk. 



Let us first comment on the behavior observed at very 
low temperature. In the quasi-classical region (A 3> T) 
all matrix elements of g tend to zero except for gVi. This 
is a general feature and is due to the fact that these 
terms are absent in the zero temperature expression. As 
previously stated the fall-off to zero is exponential, and 
in particular, for the model in exam, we have that 



g Th w T- 5 /2 e -A/T ; T<<A 
g TT « T~ 7 / 2 e- A / T . 



(18) 



Let us now look at the quantum critical region T> A, 
small temperature. The mixed term tends to a constant: 



9hT 



7T 

48 



O (T 



(19) 



6 



Instead grr must diverge at zero temperature, as it has 
to match with the diverging behavior observed in the 
ground state [B[ ■ For the diagonal term gxT one has 

9tt = ^-c v = ~+0{T), T > A. (20) 

We note in passing that this result agrees with the one for 
the specific heat obtained for general conformal theories 
[27l ] c„ = (-kcT) I (3v) , as T — > 0, since in our case the 
velocity v is one and the conformal charge c is one-half. 
We thus see that, in the present case, both ghh and <?tt 
diverge as T _1 . This has not to be the case in general, 
indeed at any quantum critical point described by a con- 
formal field theory, gxr will diverge as T _1 whereas the 
behavior of g^ is dictated by Eq. (fl~6)) . 

In this section we have analyzed the behavior of all 
the elements of the geometric tensor g. The result of this 
analysis allows one to conclude that indeed, at least for 
the specific model studied, the quantum critical and quasi 
classical regions can be clearly identified in terms of the 
markedly different temperature behavior of the geometric 
tensor g. 

B. Directions of maximal distinguishability 

The analysis carried out in the previous section can 
be further deepened by studying some useful quantities 
that can be derived by the analysis of the metric ten- 
sor g. Indeed, we will see that these quantities allow to 
give a finer description of the behavior of the system in 
the plane (h, T) and to reveal new unexpected features. 
We first start by noticing that at each point (h,T) the 
eigenvectors of the metric tensor g define the directions 
of maximal and minimal growth of the line element ds 2 B . 
Hence the vector field v*Ai(h,T) given by the eigenvec- 
tor of g related to the highest eigenvalue Am , defines at 
each point of the [h, T) plane the direction along which 
the fidelity decreases most rapidly: the latter represents 
the direction of highest distinguishability between two 
nearby Gibb's states. 

We now focus our analysis on the study of the vector 
field vm(Ii,T) in the specific case of the quantum Ising 
model, see Fig. [TJ We first observe that there clearly 
are some interesting features for small temperatures, that 
reflect the analysis previously carried out on the metric 
elements of g. On one hand, in the quasi-classical region, 
when h ~ we have that the direction of highest fidelity 
drop, is parallel to the h axis. This reflects the fact that, 
in this region, all the elements of g tend to zero except 
for the term g^. On the other hand, in the quasi critical 
region, the direction of highest fidelity drop is parallel 
to the T axis. Again, this feature can be linked to the 
previous analysis of the terms (fT5|) and f2D|) . Indeed both 
ghh an d grr diverge as T _1 , but, as (C/n 2 ) / (w/24) = 
0.70... < 1, gxT eventually becomes bigger and thus 
the direction of highest distinguishability turns parallel 
to the T axis. 











* 


■» 


2 














^ , > •* > 


v 


* 




v 




^ , ^ ^ v... "V v . 








\ 




^ -^.^^ ^\\ \ * 




v 


y 


y 


1.5 


„ „ . ^ > > > > v 


\ 


v - 


\ 


V 










\ 


\ 






V 




V 


V 


1 
















y 


v 


y 




i y y v v v v v v v v 


V 


V 


i 






H H' v 1 4 \ \\ 


w 




w 


k 


0.5 


f 1 1 \ \ \ \ t \ \ 




k 


\ 


i 

\ ' 




1 f f f f M \ \ \ \ 

1 1 1 I 1 1 \ \ \ \ > 


\ 


1 


\ 




i 1 1 i d t 1 i » 1 1 




■ 




i 





'4 \ 











0.5 1 1.5 2 
h 



Figure 1: Vector field of the eigenvector associated to the 
highest eigenvalue of g, in the plane (h, T) for the Ising model 
in transverse field. 



We proceed in our description of the phase diagram 
through the introduced vector field by examining what 
happens at the h — axis. Here one has the impression 
that a kind of singular point appear around T « 1. The 
reason for this is that at h = the system becomes the 
purely classical Ising model, which possess only classi- 
cal behavior at any temperature. This implies that the 
quantum-critical region cannot extend over this line. As 
the dispersion Afe is flat, it is straightforward to write 
down the metric tensor on the h — line. It turns 
out that g is completely diagonal meaning that eigen- 
vectors are parallel to the (h, T) axes. One sees that for 
0.852 < T < 00, g hh > g TT then for 0.101 < T < 0.852, 
9tt > 9hh and then finally, at very small temperature, 
quantum fluctuations dominates and for < T < 0.101, 
ghh > gTT- The appearance of the purely classical Ising 
line at h = 0, which forbids the quantum critical phase 
extend over this line, is related to the fact that the model 
(|14[) is invariant under the Z2 symmetry h — * —h. This 
in turns implies that the phase diagram is mirror sym- 
metric around the h = line and that there is another 
quantum critical point at h = — 1. The physical conse- 
quence is that the semiclassical ordered region is much 
smaller than one would think and the actual phase dia- 
gram is much like in Figure [5J 



Finally we now discuss another feature that can ob- 
served by studying «m(/i,T). As one can see in Fig. [2 
along the line T = h > 1 the vector field becomes par- 
allel to the vector w = (—1,1). It turns out that this 
feature can be understood analytically by studying the 
behavior of the metric tensor g when \h\ ^> 1. Indeed, 
by evaluating the dominant part of the various metric 
elements on the line T = h = t ^> 1, one sees that all 



Figure 2: Phase diagram of the Ising model in transverse 
field taking into account both critical points at h — ±1, and 
the purely classical Ising line h = 0. The arrow indicates 
the direction of highest fidelity decrease (the direction of the 
arrows is conventional, but fixed once for all) 

the Fisher-Rao terms decay as t~ 2 while ~ i~ 4 , and, 
what is most surprising, all matrix elements tend to have 
the same value in magnitude. This feature can be un- 
derstood by simply observing that when \h\ 3> 1 it is 
only the classical term proportional to the external mag- 
netic field of the Quantum Ising Hamiltonian that sur- 
vives i.e., H ~ ^Si°f ■ The density matrix of the sys- 
tem can be written as p(h,T) = exp (— hj^i °f /T)/Z; 
in this approximation the only non zero terms of the 
metric are the Fisher-Rao ones and all the covariances 
that define these terms, see ©, coincide with var(H). 
Thus, in the limit \h\ ^> 1 the Bures distance reads 
ds% = var(H)[dT 2 /T 2 - hdTdh/T 3 + h 2 dh 2 /T 4 }. If now 
one chooses the particular case T = h = t and evaluates 
the density g/L one finds that 

g(t,t) = C ( \ ~ 1 )+0(t- 3 ). 

yV ' 16cosh 2 (l/2) V - 1 1 / 

Thus, one has that on the line T — h ^> 1 the only non 
zero eigenvalue is 2var(H)/(Lt 2 ) and it corresponds to 
the eigenvector w = (— 1, 1). In this approximation, that 
amounts to neglecting the term g^ ~ i~ 4 , when moving 
along the line T = h 3> 1 i.e., along the direction defined 
by w 4 -, no changes in the state of the system occur. 



C. Crossover and metric tensor g 

We finally present some preliminary results related to 
the intriguing possibility of determining the crossover 
lines between the quasi-classical and quasi-critical region 
([T4]) through the analysis of the elements of metric tensor 
g and the induced Gaussian curvature [28j in the plane 
(h, T). The capability of the highest (in modulus) eigen- 
value of the g and of the Gaussian curvature induced by 
the metric to capture, in terms of divergencies or discon- 
tinuities, the existence of QPTs has been already tested 
in |3| and Here we would like to test whether these 
quantities are able to identify the crossover between the 
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Figure 3: Contour plot of the highest eigenvalue of g, in the 
plane (h, T) for the Ising model in transverse field. 



quasi-classical and quantum-critical region. Notice that 
the curvature of the Bures metric in the case of squeezed 
states has been studied in [2§] and an operational inter- 
pretation attempted. It is also worthwhile to stress that 
the so-called thermodynamical curvature plays a central 
role in the geometrical theory of classical phase transition 
developed by Ruppeiner and coworkers [301 ] . 

As already pointed out, at each point (h, T) the vec- 
tor field VM(h,T) defines the direction of highest dis- 
tinguishability between two nearby Gibb's states. The 
degree of distinguishability along this direction is quan- 
tified by the maximal eigenvalue Am (h,T). Since the 
quasi-classical and quantum-critical regions are charac- 
terized by significantly different physical properties, it is 
natural to investigate whether the change of the latter, 
in spite of not involving a phase transition, could be re- 
vealed by our measures of statistical distinguishability 
and by the related functionals. 

We now give a descriptive analysis of the raw data. 
In figure [3l we have plotted the contour plot of the 
maximal eigenvalue of g. The main feature is the pres- 
ence for T > of two patterns of high distinguishability 
(white) that separate the regions (h < 1,T < 0.25) and 
(h > 1,T < 0.25) from the rest of the diagram. Thus, 
the first information that can be drawn is that a change 
of parameters inside these regions implies a small change 
in the statistical properties of the corresponding ground 
states. On the contrary, if one varies h and T and moves 
from these regions towards the center of the diagram, for 
example moving along the integral lines olvM{h, T) , the 
statistical properties of the state necessarily have to sig- 
nificantly change. One can see that, the "transition" lines 
between the different regions can be extrapolated numer- 
ically by tracing the "ridge" lines of the two patterns of 
high distinguishability. It turns out that the same re- 
sult can be achieved by looking at the lines where the 
Gaussian curvature of g changes sign, see figure HI For 
example, when h > 1, one can see that along the de- 
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Figure 4: Contour plot of the Gaussian curvature of g, in 
the plane (h, T) for the Ising model in transverse field. The 
arrows indicate the zero curvature lines. 



termined transition line, T has a linear dependence on 
h - 1. 

This preliminary descriptive analysis seems thus to in- 
dicate that a neat distinction between the quasi-classical 
regions (characterized by a negative curvature) and the 
quantum-critical (characterized by a positive curvature) 
can be made on the basis of study of the metric g. This 
is indeed the first time that the use of the fidelity, and 
of the related functionals, allows to identify the crossover 
between two distinct phases. 



VI. CONCLUSIONS 

In this paper we have analyzed the relation be- 
tween quantum criticality, finite temperature and the 
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